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Abstract 

We consider the one-dimensional diffusion of a particle on a semi- 
infinite line and in a piecewise linear random potential. We first present 
a new formalism which yields an analytical expression for the Green 
function of the Fokker-Planck equation, valid for any deterministic con- 
struction of the potential profile. The force is then taken to be an asym- 
metric dichotomic process. Solving the corresponding energy-dependent 
stochastic Riccati equation in the space-asymptotic regime, we give an 
exact probabilistic description of returns to the origin. This method al- 
lows for a time-asymptotic characterization of the underlying dynamical 
phases. When the two values taken by the dichotomic force are of differ- 
ent signs, there occur trapping potential wells with a broad distribution 
of trapping times and dynamical phases may appear, depending on the 
mean force. If both values are negative, the time-asymptotic mean value 
of the probability density at the origin is proportional to the absolute 
value of the mean force. If they are both positive, traps no more exist 
and the dynamics is always normal. Problems with a shot-noise force 
and with a Gaussian white-noise force are solved as appropriate limiting 
cases. 
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1 Introduction 



The study of dynamical features of diffusion or conductivity in random en- 
vironments has been initiated in early eighties |IJ and since then intensively 
pursued @-|§. The problem is usually formulated either in a discrete form, i.e. 
by means of the Pauli Master Equation ||, assuming the transfer rates to be 
random variables, or in a continuous formulation, i.e. using the Fokker-Planck 
equation |j)-[f§ and assuming the so-called drift-function to be a stochastic 
function in space. The transport properties are obtained through an average 
over the random parameters in the equation of motion. 

The resulting dynamics in the time-asymptotic region may exhibit non- 
standard features. For instance, with a Gaussian drift function of non-zero 
mean value, a succession of dynamical phases is observed when the bias is var- 
ied. These phases are characterized by different drift and diffusion behaviours, 
normal (i.e. with a finite mobility and diffusion coefficient), or not. The ex- 
istence of anomalous dynamical phases can be traced back to the existence of 
traps with a broad distribution of trapping times ||]. 

Most of the standard treatments assume independent random transfer 
rates in the master equation |9|, 10|| , or, correspondingly, a white- noise drift 
function in the Fokker-Planck equation [fn{]-p4|. Nevertheless, as pointed out 
e.g. in [T5|], a more realistic description should take into account the possibility 
of spatial correlations in the local transport parameters flnS[|-|j20|. 

Let us for instance consider a continuous medium and imagine first that 
the particle diffuses through an array of randomly positioned minima with 
randomly distributed depths. The simplest choice for a space-correlated bias 
yielding this scenario is that of a dichotomic noise ||, with realizations al- 
ternately assuming two possible values of different signs. The potential then 
displays a succession of linear segments of random lengths and of alternately 
positive and negative slopes. One can well figure out that, in such a potential, 
the minima represent traps, some of them being deep traps with large trapping 
times. This case shares some similarities with the Gaussian-white-noise one, 
and should lead to results of the same type for the particle drift and diffusion 
properties (i.e. the existence of anomalous dynamical phases for certain val- 
ues of the parameters of the model). But, interestingly enough, the choice of 
a dichotomic noise also allows to treat the case of a quenched random force 
taking alternately two values of the same sign. Since clearly the notion of deep 
traps — and even simply of traps — makes no more sense in this case, one does 
not expect the existence of anomalous dynamical phases. 
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More specifically, in the present paper, we assume that the lengths of the 
above described constant-force segments are independent and exponentially 
distributed, in which case the quenched random force is a Markovian Poisson 
process in space. We consider a particle diffusing on a semi- infinite line, i.e. 
we impose a reflecting boundary condition at the origin. These two hypotheses 
allow for exact analytical calculation and for the discussion of a rich variety of 
physically different situations. 

The paper is organized as follows. In Section 2, our analysis begins 
with the Laplace transformation of the Fokker-Planck equation. One thus 
gets a differential equation in space, depending parametrically on the energy, 
as pictured by the Laplace variable z. The corresponding Green function in 
the presence of an arbitrary deterministic piecewise constant bias is derived. 
In Section 3, we introduce quenched disorder with general piecewise constant 
realizations. The localization probability of the particle at its (sharp) initial 
position satisfies a Riccati stochastic differential equation with a multiplicative 
noise (in space). Specifying further the noise to be the Markovian Poisson one, 
we are then able to derive a one- formula based (Eq. (|50|)) parallel analysis of the 
disorder average of the localization probability of the particle at the origin on 
the one hand and the trapping time or the time-asymptotic average velocity on 
the other hand. In Section 4, we analyze various physical situations, according 
to the sign of the mean bias and to the presence or to the absence of traps. 
Finally, Section 5 contains our conclusions. 

Generally speaking, the new results of our paper are the following. First, 
our procedure is exact for any fixed value of the Laplace variable z. On the one 
hand, this enables, at least in principle, the analysis of the disorder-averaged 
probability density at the origin for any time. On the other hand, we can 
carry out the small- z analysis and derive the exact time-asymptotic formulae 
for this quantity in various physical regimes. Another specific feature of our 
work is the application of reflecting boundary conditions at the origin. Thus, 
clearly, in contrast to the equivalent problem on an infinite line, the situation 
with a positive or a negative mean bias are not equivalent. Indeed, with a 
negative mean bias the particle is in some sense stuck to the boundary at the 
origin or pushed back towards it, while with a positive mean bias it escapes 
towards infinity, the modus of its time-asymptotic motion being controlled by 
a parameter describing the typical depth of the potential traps. 

2 Diffusion in a deterministic force 
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2.1 Homogeneous force 

Let us consider an overdamped Brownian particle acted upon by a standard 
white- noise Langevin force T(t) and by a position-dependent potential force 
Fix). Its dynamics is described by the viscous Langevin equation 



V j T x{t) = F[x(t)} + r® , (i) 



with t] being the viscosity. The correlation function of the Langevin force is 
equal to 2Dq7] 5{t — t ), where Dq = hsT/i] is the diffusion constant in the 
absence of the potential force. The corresponding Fokker-Planck equation for 
the Green function P(x, y; t) reads 



9 ~ 7s 9 



d ~ ~ Fix) ~ 
-D —P(x,y;t) + ^^P(x,y;t) 
ox rj 



(2) 



The bracketed expression represents the probability current J(x,y;t). We 
assume the initial condition P(x, y\ 0) = S(x — y) and the boundary conditions 
J(x ,y;t) = 0, J(xi,y;t) = 0. Consequently, the boundaries at x and x% are 
reflecting and the probability density is always normalized to unity. 

In order to make the following calculation more transparent, we introduce 
dimensionless variables. The potential force will be written in the form F(x) = 
Fof(x). The dimensionless coordinate is x = xFo/rjDo, and the dimensionless 
time t = £Fq/VA)- We thus get from Eq. (§): 



d d 

—P(x,y;t) = -—J(x,y;t) , (3) 
J(x,y;t) = ~P(x,y;t) + f(x)P(x,y;t). (4) 



The original density and current are connected with their dimensionless coun- 
terparts via P(x, y; t) = F P(x, y; t)/r]D and J(x, y; t) = FqJ(x, y; t)/rj 2 D. 

Let us assume for the moment a position- independent force: f(x) = f 
for x G [xo,xi]. Performing the Laplace transformation, the Fokker-Planck 
equation (H) yields a nonhomogeneous differential equation with constant co- 
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efficients, 



— -/- + * 

dx 2 dx 



P(x,y;z) = -5(x-y) . (5) 



We are using the same symbol for a given function a(t) and for its Laplace 
transform a(z) = J °° dt exp(—zt)a(t), the Laplace original (transform) be- 
ing always indicated by writing the variable t (z). Combining any particular 
solution of the nonhomogeneous equation with the general solution of the ho- 
mogeneous equation, we have P(x, y\ z) = Pn(x, y; z) + Ph(x; z), with 

P N (x,y;z) = —^—\Q(y-x)e^ a+ ^ + Q{x-y)e- {x - y)a ~ {z) \ , (6) 
2a(z) ' 

P H (x-z) = c + (z)e xa+(z) +c-(z)e- xa ~ {z) , (7) 

where a(z) = sfz + / 2 /4, a ± (z) = y/z + f 2 /4±f/2, and G(x) is the Heaviside 
function. Having acquired the general solution for the density, we calculate 
the general expression for the probability current. In the last step, the two 
functions c ± (z) are fixed from the reflect ing-boundary conditions J(x ,y; z) = 
0, J(x 1: y; z) = 0. 

The whole procedure is well known. However, in view of the following 
calculation, it is convenient to present the final result in a matrix form. We 
introduce a two-dimensional space with the basis {|1, 0) , |0, 1)} and we express 
the pair density-current as two coordinates of a single state ket: P(x,y; z) = 
(1, 0\G(x, y; z)) and J(x,y;z) = (0, l\G(x, y; z)). Adopting this convention, 
the result of the present simple example reads 

\G(x, y; z)) = W(xi - x; z)\l, 0) T(y; z) - Q(y - x)\N(y - x; z)\0, 1) . (8) 

Here we have introduced the abbreviation 



L[V ' Z) (0,l\yN( Xl -x ;z)\l,0)' {) 
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and the matrix 



\N(x;z) 



I a -( z )e xa ~^ + g + (z)e- xa+{z) e xa ~ (z) - e - XQ+{z) \ 

2a(z) 2a(z) 
e xa~(z) _ e -xa+(z) a + {z)e xa ~ {z) + oT (z)e~ xa+ {z) 



\ Z 2a(z) 2a(z) J 

(10) 

Notice that, at given z, W(x; z) satisfies the "dynamical equation" 



^-\N(x;z) = H(z)\N(x;z) , H(z) = [ ^ * ) , (11) 




with the position x playing here the role of time. Eq. (p~Q|) gives W(0; z) = I 
(unity matrix), i.e. we have formally W(x; z) = exp[xH(z)}. 

Eq. (|j) yields the complete picture of the resulting motion (i.e. we can 
compute P(x, y; z), J(x, y; z) and these functions can be inverted into the time 



domain [23]). For instance, taking x = 0, x\ = I, y — > + , and x = 0, the 



probability density at the origin emerges as the ratio of two matrix elements: 



Pm n r\ - Mm\z)\l,0) _ 1 a-(*)e to -('> + a + (z)e- to+ <«> 
K ' ' 1 (0,l|W(Z;z)|l,0) z e ' a "W-e- fa+ W ' 



Moreover, for the semi-infinite line, we have lim^oo P(0, 0; z) = a~{z)/z. In 
this case, we get following picture. Having / < 0, the force pushes the diffusing 
particle against the boundary. In this case, the time-asymptotic value of the 
probability density at the origin is |/|, the asymptotic value of the mean par- 
ticle position is i.e. the time-asymptotic velocity is zero. On the other 
hand, when / > 0, lim^oo P(0, 0; t) decreases exponentially to zero and the 
mean position increases linearly with time, the velocity being just /. Finally, 
in the marginal case / = 0, P(0, 0;£) behaves asymptotically as 1/y/nt, the 
mean position increases as ^/ttI, and the asymptotic velocity is zero. 



2.2 Piecewise constant force 

Let the original interval [xo,xat] be divided into iV segments [xk-i,Xk], k = 
1, . . . , N, with lengths If. = x^ — Xk-i- Let be the constant force in the k-th 
subinterval. We assume that the particle has been originally placed in the 
M-th segment, i.e. P(x,y;0) = 8(x — y) with y G [xm-i,^m]- The boundary 
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conditions at xo and xn are again reflecting. 

The procedure for solving the Fokker-Planck equation will be parallel 
to that in the above simple example. The general solution in the M-th seg- 
ment consists of two parts. First, the particular solution of the nonhomo- 
geneous equation (|) (with f M instead of /) will assume the form @ with 
olm{z) — \J z + /m /4 instead of a(z) and afj(z) = olm{z) ± /m/2 instead of 
a (z). Second, the general solution of the homogeneous equation in the M- 
th subinterval assumes the form (|7|), again with (^(z) instead of a ± (-z), and 
with two arbitrary functions c^-(z) instead of (^(z). The general solution in 
the A;-th subinterval, k ^ M, is also of the form (^) with the substitutions 
a ± (z) — > af(z) and cr^z) — > c k (z). Altogether, the whole general solution 
depends on the 2N functions c^{z), k = 1, . . . , N. These are fixed from the re- 
quirements J(x ,y; z) = J(x]y,y] z) = at the reflecting boundaries and from 
the continuity conditions for the probability density and for the probability 
current at the intermediate points x±, X2, ■ ■ ■ , Xjv-i- 

The final result can be again expressed in matrix form. We designate 
the "evolution operator" for the k-th segment as \Nk(x;z) — it is defined by 
the expression (|IUD with the substitutions a(z) — > aik(z) and a ± (z) — > cn^(z). 
Further, we introduce the notations 



— W m (/ m ; z)\N m+ x{l m+ i] z) . . . \N n (l n ; z) , (13) 

r , v (0, l|Wi, M -lW M (y - X M -l] 2T)|0, 1) 

TmAv]z) = (MIW^O) • (14) 



Finally, let \Gk{x, y\ z)) be the value of the state ket in the k-th segment, that 
is \G(x,y;z)) = \Gk(x,y; z)) for x G [xk-i,Xk\- The final result of the present 
Section then reads 

\G x {x,y;z)) = Wi(xi - x; z)W 2 ,at|1, 0) T M>N (y; z) - 

- Wi(xi -x;z)\N 2M -i^M(y-x M ~i;z)\0,l) , (15) 

\G M -i(x,y; z)) = \N M -i(x M -i - x; z)\N M<N \l, 0) T M)N (y] z) - 

- W M _i(x M -i - x; z)\N M (y - x M -i, z)\0, 1) , (16) 
\G M {x,y,z)) = \N m {xm - x; z)\N M+liN \l,0) T M , N (y; z) - 
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- Q(y-x)\N M (y-x;z)\0,l) , (17) 
\G M +i(x,y;z)) = \N M+ i(x M +i - x; z)\N m +2,n\ 1 i °) r M,iv(y; z) , (18) 



\G N (x,y;z)) 



\N N (x N -x;z)\l,0)T MjN (y;z) . 



(19) 



It is easy to check that the boundary conditions and the continuity conditions 
are actually satisfied and that the resulting probability density is properly 
normalized. 

Having at hand the complete Green function for the given composition 
of the segments, we now proceed to the analysis of some consequences. First, 
let x = 0, xjv = I, y + , and x — 0. The probability density at the origin 
assumes a particularly simple form: 



P(0,0;z) = 



(l,0|W 1 (/ 1 ^)W 2 (/ 2 ;z)...W jV (/ jV ^)|l,0) 
(0,l|W 1 (Z 1 ;z)W 2 (Z 2 ;z)...W JV (Z JV ;z)|l,0) ' 



(20) 



In Section 3, this expression will be taken as a starting point for the disordered- 
medium calculation. 

Second, consider the value of the probability density at point y in the 
M-th segment, i.e. at the initial position of the particle. Taking xo = and 
x — > y + , P(x, x; z) equals the ratio 



(0, 1|Wi, M -iW m Qe - xm-iI z)\0, 1) (1, 0|W m (^m - x; z)W M +i,jy|l, 0) 
(0, ljW^Zi; z)... W m (/m; z)... \N N (l N ; z)\l, 0) 



(21) 



Due to the exponential nature of the operator \Nm(Im] z) in the denominator, 
we can split it as a product of two factors, Wm(x — xm-i] z)\N m(xm — x;z). 
Further, we can insert in between the resolution of the unity operator, which 
yields 



P(x,x;z) = 



+ 



(0, l\\N k (x k - x- z)\N k+1 (l k+1 ; z) . ..\N N (l N ; z)\l, 0) 
(1, 0\W k (x k - x- z)W k+1 (l k+1 ; z)... \N N (l N ; z)\l, 0) 

(0, l|Wi(Zi; z)... Wfc-^Zfc-i; z)\N k (x - x k . i; z)\l, 0)" 
(0, ljW^Zn z)... \N k _ 1 (l k _ 1 ; z)\N k (x - x k _ i; z)\0, 1)_ 



(22) 
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For any finite x, the small- ,2 limit of the second term turns out to be zero. Thus 
the quantity lim 2 ^ + P( x , x \ z ) is equal to the small- z limit of the density at 
origin for a new medium of the total length / — x. The new interval consists of 
N— k+1 segments of the lengths x k —x, l k+ i, ■ ■ ■ , In, the constant forces in these 
segments being fk, ■ ■ ■ , In- This simple result will be of considerable value in 
the analysis of the disordered medium. Namely, imagine that the particle is 
launched from the point x. The total time T(x) spent in the interval (x, x+dx) 
is equal to 



T(x) dx = 



/ P(x,x;t)dt 
Jo 



dx = lim P(x, x; z) dx , (23) 

z— >0+ 



i.e. T(x) can be related to the probability density at the origin for the above 
mentioned new interval. Of course, T(x) can only be finite for an interval of 
infinite length. However, even for the semi-infinite interval, T(x) is finite only 
if the particle can escape to infinity, that is, if the force in the last segment 
(which is necessarily of infinite length) is non-negative. 

Third, let us consider the thermally-averaged position of the particle 
M(l; t) (we take again x = 0, xn = I, and y — > + ). Introducing the Laplace 
transform M(l; z) = Jq dxxP(x, 0; z) and integrating the Laplace transform of 
the Fokker-Planck equation, we have 



i r i 

M(l;z 



x [ X dx'P{x', 0;z) - f dx [* dx' P{x', 0; z) = 
Jo Jo io Jo 

I r l 1 1 r l 

= - - / dx-\l - J(x,0;z)] = - / J(x, 0; z) dx . (24) 
z Jo z z Jo 



Thereupon, taking the projection of the above Green function, we can write 



* rx k 

(0,l\\N k (x k -x;z)\N k+1 (l k+1 ;z)...\N N (l N ;z)\0,l)dx 

{ ' ,Z) = ~z (0,l|W 1 (/ 1 ^)W 2 (/ 2 ^)...W JV (/ J v^)|0,l) ' 

(25) 

In the next Section, this result will allow to connect the thermally-averaged 
position M(l; z) with the Laplace transform P(0, 0; z) of the probability density 
at the origin. We now turn to the detailed analysis of this latter quantity in 
the disordered medium. 
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3 Piecewise constant random force 



The preceding calculation in valid for any deterministic piecewise constant 
force. Such a force can be conceived as a member of a randomly constructed 
family of functions, that IS, clS 8b realization of a stochastic process. An arbitrary 
fixed realization is associated with a given weight. Consequently, the same 
weight is attributed to the probability density at the origin P(0, 0; z) for this 
realization. We are thus guided to the question: what is the probability density 
of the random variable P(0, 0;z)? This Section presents the exact answer 
for a semi-infinite medium and a special type of the stochastic process to be 
fully specified in Subsection 3.3. In the first and the second Subsection, our 
reasoning is valid for any piecewise constant random force. 

3.1 Stochastic Riccati equation 

Let us consider the following system of two stochastic differential equations, 

^(X; z)} = H(X; z)m\; z)} , H(A;z)=("^ jj, (26) 

where 0(A), A > 0, is a piecewise constant random process. Let us introduce 
the projections R(X; z) = (l,0\ip(X; z)), S(X;z) = (0,l\ip(X; z)), and take the 
initial conditions R(0; z) = 1, S(0; z) = 0. The system (^6|) can be solved for 
any specific realization of 0(A). Actually, consider the composition described 
at the beginning of Subsection 2.2. (c.f. also Fig. 1), and let the "time" A 
equal / — x. We have 0(A; z) = fx for A G [0, In]- If A increases, the evolution 
is controlled by the operator \Nn{X; z). At the end of this interval, the state 
\N n (In', ^)|1)0) represents the initial condition for the semigroup evolution in 
the succeeding interval A G [ijv, ijv-i + In]- This evolution is governed by the 
operator \Nn-i(X; z). Repeating this reasoning, the solution at the end of the 
iV-th interval, i.e. at the point A = /, reads 

R{l-z) = (l,0|W 1 (/ 1 ^)W 2 (/ 2 ;z)...W^(/ iV ^)|l,0) , (27) 
S(l-z) = (0,l|W 1 (Z 1 ;^)W 2 (Z 2 ;z)...W JV (Z J v; 2 ;)|l,0) . (28) 

By comparing with Eq. (PU[), it can be seen that the random variable P(0, 0; z) 
is identical to the ratio R(l; z)/S(l; z). From now on, the random variable 
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P(0, 0; z) will be designated as P(X;z). We introduce its probability density 
7r(p, A; z) = ( 5[p — P(X; z)] ), where the symbol (. . .) denotes the average over 
the quenched disorder. Taking the projections of Eq. (p6|), the random vari- 
ables R(X; z) and ^(A; z) obey the system of stochastic differential equations 



-^R(X; z) = -0(A) J2(A; z) + S(X; z) , ^-S(X; z) = zR(X; z) , (29) 



with A representing the total length of the interval accessible to the diffusing 
particle. Finally, on comparing the A-derivative of the product R(X; z) = 



P(X; z) S(X; z) with the first equation ( p9|) and on dividing by S(X;z), we 



obtain the stochastic Riccati differential equation obeyed by P(X; z): 

^-P(X-z) = -zP 2 (X;z)-<P(X)P(X-z) + l , P(0;z) = +oo. (30) 
dX 

For purely operational reasons, we introduce the random variable Q(X; z) = 
zP(X; z). One gets from Eq. (|30|): 



4rQ(X;z) = -Q 2 (X;z)-(f>(X)Q(X;z) + z , Q(0; z) = +oo . (31) 
dX 



The density n(q, A; z) = ( 5[q — Q(X] z)] ) is related with the density ir(p, A; z) 
by 7r(p, A; z) = zn(zp, A; z). Regarding the relation 



limP(0,0;t) = lim zP(X; z) = lim Q(X; z) , (32) 



the small-2; limit of the density n(q, A; z) describes the time-asymptotic prop- 
erties of the random variable P(0, 0; t). 

3.2 Mean trapping time and mean velocity 



Let us now return to the reasoning associated with Eqs. (pil)-(^). We shall 
call the variable T(x) = lim z ^ + P(x,x; z) the trapping timef} Presently, the 
density P(x, x; z) and hence also the trapping time are random variables. Per- 

§We shall use this term even if the potential is monotonous, although a more appropriate 
name in this case would be the dwelling time. 
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forming the small- z limit in Eq. (|22|) , we have connected the trapping time with 
the particle-position probability density at the origin of an interval shorter than 
the original one. However, if the original interval is of infinite length, the same 
is true for the new one. Since the quenched random force is described by a 
stationary process, the probability density at the beginning of the new interval 
is stochastically equivalent with the density at the beginning of the original 
one. Summing up, one has: 



T(x) = lim lim P(x, x; z) = lim lim P(X; z) . (33) 

z^0+ A^oo z^0+ A^oo 



The trapping time is position-independent, i.e. T(x) = T, its mean value 
being 



r == lim lim(P(A;;z)) = lim lim -(Q(X;z)) . 

2^0+ A^oo ' 2—>0+ A-^oo z 



(34) 



In the preceding Section, we have introduced the Laplace transform of 
the thermally-averaged position, M(X;z). Presently, it is again a random 
variable. Using the definition of S(X;z), the probability current in Eq. 
can be rewritten as J(x, 0; z) = S(X — x; z)/S(X; z), which yields 



M{X;z) 



I x | / dxS(x;z) 

- / a!xJ(x,0 ] z) = - Jo 

z Jo z o(\;z) 



(35) 



where we have used the stationarity of the quenched random force. The deriva- 
tive of this equation yields 



-^-M(A; z) = -- zP(X; z)M{\; z) 
d\ z 



M(0;z)=0. (36) 



Finally, introducing the Laplace transform of the thermally-averaged velocity, 
V(X;z) = zM(X;z), the time-asymptotic velocity for the semi-infinite line is 
given by 



lim lim V(X;t) = lim lim z / dA'exp 

t— >oo A— >oo 2— >0+ A-+oo JQ 







/ dX' exp 


-z [ X dX"P(X"; z) 


/o 


J\> 



(37) 
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where we have used lim^oo lirriA-Kxj ^(A; t) = lim 2: _ i .o+ I^a-kx) zV(X; z). The 



last formula can be rewritten in the form which reveals the well known |4], [12 
self- averaging property of the time-asymptotic velocity, when it is nonzero. 
Indeed, assuming that the limit lim 2 ^ + hm^_ >0O (P(A; z)) is finite, we can write 



lim lim V(A; t) = lim lim z \ dX' exp \—z(X - A )(P(A; z))] x 

t— >oo A— >oo z^0+ A^oo JO 



x expi-z j\x!'[P{X";z) - (P(A;z))]j 



(38) 



Due to the above assumption, if A —>■ oo, the integral in the second exponent 
represents a random variable which is finite (for high enough A", the typical 
trajectory of P(X";z) swings around the mean value (P(X;z))). Thereupon, 
for small z, the second exponential tends to unity and the remaining integra- 
tion yields the reciprocal value of the non-random number fl34|) . Thus, when 
the mean trapping time r is finite, the asymptotic velocity is a self-averaging 
quantity equal to r _1 p[ If the mean trapping time diverges, the asymptotic 
velocity vanishes. In this latter case, the disorder-averaged time-asymptotic 
mean position either tends to a constant (for a negative mean force), or in- 
creases slower than linearly. 

For the sake of completeness, the second thermally-averaged moment 
z) can be also connected to the probability density 
at the origin. Actually, first, the Fokker-Planck equation implies iV(A; z) = 
(2/z) / A dxxJ(x, 0; z). Thereupon, on deriving this expression, we get 



-jrN(X; z) = 2M (A; z) - zP(X; z)N(X; z) , iV(0; z) = . (39) 
d\ 



Summing up, the first and the second thermally-averaged moments obey a 
system of stochastic differential equations (p6|), (p9|), with P(X; z) playing the 
role of the "input" noise. 



^For the exact proof of this statement one needs the simultaneous probability density 
i/j(p,w, A; z) for the random variables P(\;z) and W(X;z) — z 2 M(X;z). Assuming the 
condition in the text, one has to prove lim z ^ + li m A— »oo / V'CPj w > z ) dp — S(w — t~ x ). 
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3.3 Dichotomic random force 

Let the forces in the individual segments assume alternately just two values, 
F± = Fof±. We shall always take /_ < /+; the equality would imply a 
position-independent constant force. Let the lengths of the constant-force 
segments be independent random variables. The generic probability density 
for the (dimensionless) lengths of the constant-force segments will be taken 
of the form p±(A) = n± exp(— \n±), where l/n± denotes the mean length of 
the segments with the force f±. Due to this assumption, the resulting four- 
parameter stochastic process <f)(\) is Markovian and it usually referred to as 
the asymmetric dichotomic noise [§]. We shall always work with a stationary 
dichotomic noise. We can set 

pg W A)) = / -** + / +"- , W (A)0(A')> = fexp(-I^T) , (40) 

where we have introduced the intensity a = f n_n + (f + — /-) 2 / (n_ + n+) 3 , and 
the correlation length A c = (n_ + n+) -1 . The statistical properties of the 
stationary noise are invariant with respect to the inversion A — ► —A and to the 
translation A — > A — I. 

Let us now focus on Eq. ([31]). It is convenient to associate with the 
variable Q(X; z) an overdamped motion of a hypothetical particle. While the 
"time" A increases, this particle moves alternately under the influence of the 
"forces" 

K±(q; z) = - q 2 -f ±q + z = -[q- q ± {z)} [q - q' ± (z)] , (41) 

where we have introduced the four quantities 

q±(z) = ^z + p ± /4-f ± /2 , q f ± (z) = -y/z + f£/4-f ± /2. (42) 

Notice the ordering q' + (z) < q'_(z) < < q+(z) < q~(z), valid for any real 
positive z and for any values of the parameters /_</+. The corresponding 
"potentials" U±(q; z) = — J K±(q; z) dq display minima at q±(z) and maxima 
at q'±(z). Starting from its initial "position" at infinity, the particle always 
slides either towards q~(z) or towards q+(z) < Q-(z). For any fixed A the par- 
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ticle can only be found between the coordinate valid for the potential U + (q; z) 
and that valid for U-(q; z). Accordingly, for an arbitrary nonzero z, the prob- 
ability density n(q, A; z) vanishes outside the finite interval 

- q + (z)e x ^-q' + (z)e x ^ q.(z)e x ^ - g'_{z)e x ^ ] 

e A?+W _ e V+W ' e \q-(z) _ e Xq'_(z) J ' ^ 

and it displays two 5-function contributions at the edges of this support, their 
weights being n T exp(— An±)/(n_ + n+). This singular part describes an ex- 
ponentially decreasing probability of having just one segment in the whole 
interval of length A. Obviously, in the limit A — ► oo, the support is simply 
[q+(z), q~(z)] and the singular part is missing. 

In order to solve Eq. ( |31| ) for the dichotomic noise in question we follow 
the standard steps as described in ||. First, we introduce the joint densities 



K±(q, A; z) dq = Prob { Q(X; z) e (q,q + dq) and 0(A) = / ± } . (44) 



One has n(q, A; z) = «_(g, A; z) + K+(q, A; z). These densities obey the coupled 
partial differential equations 



d 






d 


dX 


_ K + (q,\;z) _ 




dq 
















— n_ 



K-(q;z)K-(q,\;z) 
K + (q;z)K + (q,X;z) 



n i 



K+(q,X;z) 



(45) 



We are looking for the stationary solution K±(q',z) = 
Hence we remove the A-derivative on the l.h.s. of Eq. 
two new functions 



limA^oo K±(q, A; z)Q. 
Introducing the 



= r^K + ( q ;z) + nX{q; z) Z) " n+Mg; * )] ' (46) 

rj(q;z) = K_(q; z)K-(q; z) + K + (q; z)n + (q; z) , (47) 

**The nonstationary case would describe the diffusion on a finite interval, its analysis 
being very difficult even in a simpler context |2^| and of minor physical importance. 
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and carrying out the corresponding substitution in Eq. (|45|), we arrive at two 
independent equations: 



1 d 

£(q; z) dq 



£(?;*) 



n. 



K-{q;z) K+(q;z) 



(4* 



Hence the function r](q; z) is simply a constant, which in fact is equal to zero. 
Actually, Eq. (|4T| ) yields K±[q±(z); z] = 0, and conservation of probability 
entails K±[q T (z); z] = 0. As for the differential equation for £(q; z), its solution 
reads 



C(z) 



q-(z) - q 

q-qL(z) 



"-(*) 



q - q+{z) 



q-q+(z) 



u + (z) 



V±iZ 



n± 



Az + p ± 



(49) 

where C(z) is a normalization constant. Finally, inverting the transformation 
{([HP, (f47j)}, we get K±(q;z) = ^(,(q; z) / K±(q; z) and the stationary density 
K,(q; z) = limA^oo A; z) reads 



K(q;z) 



1 



+ 



1 



x 



C(z){[q.(z)-q][q-qL(z)} [q - q + (z)\ [q - q' + (z)\ 

u+(z) 

&[q;q+(z),q-(z)] > 



X 











q-(z) - q 


q 


-q+( z ) 


q - q'_(z)_ 




q 


-q' + (z)\ 



(50) 



where we have denoted^ 0(g; x, y) = Q(q — x)Q(y — q). In the final step, the 
const ar 
that is 



constant C(z) is fixed from the normalization condition fjj $ n(q~, z) dq = 1, 



C(z) = (/+-/_) 



-W [q-(z) - qf-W" 1 [q - 



I q+{Z ) q [q-q'-(z)]»-W+ 1 [q-q' + 

Some steps of the underlying calculation are given in the Appendix — the re- 

^The product of the two Heaviside functions guarantees the vanishing of the probability 
density outside the interval [q+(z), 



dq ■ (51) 
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suiting form of the above integral can be written as 



= 1 [q-{z) - q+(z)]"-^y+^ 



.J4z + 2 (A)) [?+(*) ~ q'-{z)Y-W[q_{z) - <f+(z)]"+W 
x B[v-(z), v-(z)\ T v+(z), v-(z) + v+(z) + 1; . (52) 

Here £>(x, j/) = T(x)T(y)/T(x + y) denotes the Euler beta function, J 7 (a, b, c; x) 
is the Gauss hypergeometric function [F£J, and we have used the abbreviation 



..( 7 \ = - q + (z)][q' + {z) - q' + (z)} 



Eqs. (0), ( p2|) for the stationary density represent the main result of the 
present Section. The corresponding moments can be obtained by the usual 
integration. Again, the ensuing calculation is commented on in the Appendix. 

4 Discussion 

Our general four-parameter description of the dichotomic force provides a rich 
spectrum of special regimes, which can be analyzed using Eq. fl50D, where, as 
indicated above, the probability density n{q\ z) stands for liniA^oo n(q, A; z). In 
the same way, we shall use the simpler designations Q(z), n(p; z) and P(z) for 
the stationary values Q(X; z), 7r(p, A; z) and P(X; z). 

4.1 Both forces are negative (/ _ < / + < 0) 



The mean force \i in Eq. fl40|) is negative and the potential consists of segments 
with a positive slope (c.f. Fig. 2). For any arbitrary realization of the quenched 
noise, the particle cannot escape to infinity — it can be found with probability 
one in a finite region. Intuitively, one expects a nonzero time-asymptotic mean 
value of the probability density at the origin, and a finite time-asymptotic value 
of the thermally-averaged mean position. 

Consider the small-z limit of the probability density fl5"0|) . First, one has 
lim 2 ^ + u±(z) = n±/\f±\. Second, the small-z limits of the expressions q±(z) 
are |/±|, whereas q±(z) behave as —z/\f±\. Analyzing the expressions in (|43|), 
the support of the probability density lim 2 ^ + n{q\ z) is the interval [|/+|, |/-|]- 
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Finally, we have u(z) — > 0, i.e. the hypergeometric function in Eq. ( j5^) tends 
to unity. On collecting these observations, one gets 



lim /c(p;z) = \n\ ^ 1 + ^nJ » # 1 ( Tf~ v TJ 



(IJ-l-^^-li+l)^" 1 



ri_ n i 



e(g;|/ + |,|/_|) . (54) 



The corresponding moments lim z ^ + \Q k ( z )/ are & U finite and they can be 
computed analytically by direct integration. In particular for k — 1 one gets 
lim 2 ^ + (Q( z )) — M- Thus the time-asymptotic mean value lim t ^ oo (P(0, 0; t)) 
in the semi-infinite line is seen to be equal to the absolute value of the mean 
force, \fi\. Note that this result cannot be obtained from the solution of the 
corresponding free-diffusion model on an infinite line (i.e. without the reflect- 
ing boundary condition at the origin). The higher moments are not so simply 
related to the properties of the random force. Further, in the present case, the 



mean trapping time (|3H ) is infinite and the time-asymptotic velocity vanishes. 
The particle is in some sense stuck to the origin. If /_ < / + = 0, a slightly 
more complicated calculation reveals the same general conclusions. 

4.2 Forces are of different signs (/ _ < < /+) 

The potential forms a system of traps (c.f. Fig. 1). The traps can only be effi- 
cient if the ratio n + /f + is comparable with the ratio n_/|/_|. Otherwise, they 
are typically "shallow" and they do not represent sufficiently effective obsta- 
cles for the particle motion. The "trap-permeability" parameter, as defined by 
9 = f — n + /f +1 will play an important role in the following discussion 



PJ, [TjJ. In fact, it is proportional to the mean force: 9 = /^(n_ + n+)/|/_|/+. 
Having fixed the forces f±, the mean force \x can be either positive or negative, 
depending on the parameters n±, the value /i = separating two regions with 
essentially different time-asymptotic dynamics. 

Let us consider again the small- z limit of Eq. (|50| ). The quantities q~(z) 
and q'_{z) behave as in the preceding Subsection. Presently, however, one has 
q+(z) ~ z/f+ and q' + (z) — > — /+. Thereupon, the small- z limit of the support 
3D is now the interval [0, |/_|]. Further, the variable fl53D diverges and one 



has to use the analytic continuation of the hypergeometric function in Eq. (|5^). 
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4.2.1 Negative mean force (// < 0) 

In the small- z limit the normalization constant C(z) alone converges to a finite 
number and we can safely carry out this limit separately in C(z) and in the 
rest of the expression (p0|). The result reads 



\f-V+ B 



1 










'/+ 





x g £ ^ 'Qf^-q)^ '(q + U) 1 Q{ q] ^\f_\) . (55) 



All the moments of this limiting density exist and can be computed analyti- 
cally. Let us just quote here the result for k — 1: lim 2 ^ + (Q( z )) — | A* | , i-e. we 
have again lim^ oo (P(0, 0; £)) = The mean trapping time diverges and the 
time-asymptotic velocity vanishes. As compared to the previous Subsection, 
the presence of traps does not modify the modus of the asymptotic dynamics. 

4.2.2 Zero mean force (// = 0) 

In this Sinai-like case, the small- z analysis of the general expression for the 
first moment (P{z)) together with the Tauber theorem for the inverse Laplace 
transformation |26| yield the logarithmic decay 



^(0,0;;))'-°°^-^--^. (56) 
x v " n_ +n + log* 



4.2.3 Positive mean force (// > 0) 

In the small-z limit, the normalization constant C(z) diverges as z~ e , 6 > 
0. The limiting density flST^) is concentrated at one point: lim 2 ^ + n(q] z) = 
S(q). All the moments Q k {z) are self-averaging quantities, their (non-random) 
limiting value being zero. On the other hand, the density lim z _ > o+ ?r(p; z) is 
a well behaved function, concentrated on the interval +oo[. Actually, if 
we first handle the substitution n(p; z) = ZK,(zp; z) in Eq. ( [50"D and afterwards 
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exercise the small- z limit, we get 



n_ n\ 



hmirfaz) = — — B [ m ~ ~T I x 

n_+n + - . tttt-^ \f+ \f-\ f+ J 

x " ( p "£)' + e K' + °°) (57) 



Here we come to an essential conclusion |fL9fl : the moment lim 2 _ >0 + (P k { z )) 
is only finite if fc < 6 1 . Specifically, for 9 e]0, 1[, the limit lim z _ >0 + (P(z)) is 
infinite, i.e. the mean trapping time is also infinite and the time-asymptotic 
velocity is zero. More precisely, using again the Tauber theorem, we get 



,pm 0-t)) 'IS? == r Hl/-l) r\M±Ag I (58) 

{F{W >> ~ „ + („_+„ + ) r2 ^ rw i/.r f (58) 

If # > 1, the first moment lim z ^ + (P{ z )) is finite and its reciprocal gives the 
(self-averaging) time-asymptotic velocity: 



lim lim V(A;f) = (0-l) 7 (59) 

i-*ooA->oo v 1 v y (n_ + n + ) 2 - n_|/_| + n + f + v 7 



However, if 6 e [1, 2], the small- z limit of the second moment lim 2 _ >0 + (P 2 (z)) 
is infinite. This can be shown to imply the vanishing of the (static) diffusion 
constant for the disorder averaged dynamics (this quantity is not discussed 
here) . 

4.3 Both forces are positive (0 < / < /+) 

In this case, the slope of the potential is always negative: the particle just slides 
towards infinity. The small- z limit of the stationary probability density for the 



random variable P(z) again follows from Eqs. (50), (p2|). The normalization 



constant alone tends to zero and one must operate with the whole expression 
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(|50D . The result reads 



lim 7r (p; z) 

z— >0+ 




i> -^-\ + e( P ;^,j-) . (60) 



/+/ V7+7- 

Obviously enough, all the moments of this density exist. For k = 1, we get 



t = hm (P (z)) = v — . 61 



As expected, the mean trapping time is finite. The time-asymptotic velocity is 



self- averaging and equals the reciprocal of the above expression |L8], |19[. There 
is no anomalous dynamical phase in this case. 

4.4 White shot-noise limit 

Originally, the dichotomic quenched force has been described by four param- 
eters, n± > and f±. Another convenient equivalent four-parameter set is 
the mean force /x, the intensity a, the correlation length A c = l/(n_ + n+), 



already introduced in Eq. (f40|), and the "non-Gaussianity" parameter |2| 



7 = |/+ — /_|/(n_ + n + ), the meaning of which being explained below. It 



is well known H, 21 that an appropriate limit of the dichotomic noise yields 



the Poisson white shot-noise. Actually, consider the parametrization 



/- = > /+ = i > n - = -2 > n + = - ■ ( 62 ) 



If we increase the parameter £, the force / + increases and the mean length 
of the segments with the force /+ tends to zero. In the limit £ — > oo, the 
parameters /i, a, and 7 keep their values, whereas the correlation length A c 
tends to zero. The limiting form of the correlation function in Eq. ( flUl) is 
2<T 5(X — A'). After the indicated limit, the quenched force displays an array 
of randomly positioned 5-impulses on the constant background — (a — 7/x)/7 
(c.f. Fig. 3). The mean (dimensionless) distance between the (^-impulses is 
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7 2 /cr, their weights being randomly distributed with the probability density 
7 _1 exp(— w/j)Q(w). Thus the parameter 7 represents the mean weight of the 
impulses. On the whole, in the present Subsection, the random potential is 
described by the three parameters a > 0, 7 > 0, and fi. The potential wells 
only exist if /_ < 0, i.e. if 7// < a; Fig. 3 illustrates the typical form of the 
potential in this case. 

Let us now carry out this limiting process in Eqs. (|50[) and fl52|). We get 
q+(z) — > and u+(z) — > I/7, i.e. the support of the density K(q;z) comes to 
be the interval [0, q~(z)]. The density itself reads 



K(q; z) 



x T~ 



v-(z), v-(z) + 1, v-{z) + 



7' g_(z) - gi(z) 



x 



x q lh[ q _( z ) - q Y-W-i[ q - q i( z )]-»- W-i 6 [g; 0, ?_(«)] . (63) 



We are again interested in the small-z limit of this probability density and 
in its moments. We shall restrict the discussion to the physically interest- 
ing case with traps, i.e. 7/z < a. In this case, the variable of the Gauss 
hypergeometric function in Eq. (^) tends to 1~ and we use an appropriate 
analytic-continuation formula [^3[ . 



If fi < 0, the time-asymptotic value of the averaged density at origin is 
again lim^oo (P(0, 0;£)) = If \i = is zero, we observe the logarithmic 
decay 

(P^O-^^a^. (64) 



If fi > 0, the time-asymptotics is controlled by the trap-permeability parameter 
9. Presently, it can be written as 9 = — 7x7). For 9 e]0, 1], the trapping 
time diverges and the asymptotic velocity vanishes. More precisely, we have 
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Finally, if 9 > 1, the first moment lim z ^ + (P(z)) is finite and its reciprocal 
gives the (self-averaging) time-asymptotic velocity: 



lim lim V(X;t) = (6-1 

t—xx A— »oo 



1 + 7 



(66) 



4.5 Gaussian-white- noise limit 



It is well known [pi |2T| that an appropriate limit of the dichotomic noise yields 



the Gaussian white noise. However, the Gaussian white noise can be also 
obtained as a limit of the above- described white shot noise: one simply sets 
7 — ► + . This means that the mean weight of the 5- impulses of the force 
tends to zero, and simultaneously their density n_ = a/7 2 increases, such that 
the product (density) x (mean weight) 2 remains constant. The bias \x and the 
intensity a keep their values and the correlation function in Eq. ( fHf ) is again 
2<7<5(A — A'). The quenched force displays an infinitely dense array of S peaks 
in both directions, their weights being infinitely small. Notice that this limit 
can only be achieved if we start with the dichotomic random force taking two 
values of different signs. 

The Gaussian-white-noise results simply emerge after we carry out the 
small-7 limit in Eqs. fl6"3"|). Particularly, we have q~(z) — > 00, i.e. the support 
of the probability density n(q; z) becomes the infinite interval [0, +00 [. The 
density itself reads 



K(q; z) 



z 



6/2 



1Q, [ - 

a 



T exp 




9(g;0,+oo) 



(67) 



in accordance with the result found in Refs. j|, |T!J. Presently, the trap- 
permeability parameter simply measures the ratio between the mean force 
and its intensity: 9 = fi/a. 

The formula ( |6^) is valid for arbitrary values of the parameters a and 
\i. If fj, < 0, we have again lim^oo (P(0, 0;t)) = In the Sinai case, i.e. for 
H = 0, the asymptotic behaviour is again given by Eq. (|64|). The corresponding 
result for the infinite line without the reflecting boundary condition at the 
origin is (P(0,0;t)) w (r/(\ogt) 2 . Thus the presence of the boundary 
slows down the decay of the disorder-averaged probability density at the origin. 
Having fi > and 9 g]0, 1[, the asymptotic velocity vanishes and the averaged 
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probability density at the origin decreases algebraically as 



For example, # = 1/2 yields the exact solution (P(0,0;t)) = l/v7rt, valid 
for any time. Finally, when > 1, the time-asymptotic disorder- averaged 
mean position linearly increases, the self- averaging velocity being (8 — l)a. 
The damping of the disorder-averaged probability density at the origin can be 
exemplified by taking 6 = 3/2: we then get (P(z)) = 2/(cr + that is 

< p(0 - 0i «» = h - \ exp (r H ) erfc (H ^ ^ ^ ■ (69) 



This asymptotic behaviour should be contrasted with the exponential damping 
which takes place in the presence of a positive homogeneous deterministic force, 
as found in Subsection 2.1. 

5 Conclusion 

In the present paper, a transfer-matrix-like method for solving diffusion prob- 
lems in a piecewise linear random potential has been introduced. The formulae 
for the Green function derived in the second Section can be easily adapted to 
numerical simulation of the diffusive motion in any potential of the mentioned 
type. For example, the force can be assumed to be a semi-Markovian or a 



non-Markovian variant of the dichotomic noise ||27||-||28||, it can exhibit jumps 



of random magnitudes (kangaroo process [26]), etc. For any such process, our 
analysis is valid up to Subsection 3.3. Our subsequent choice of a Markovian 
dichotomic process has been dictated by a relatively direct possibility to get 
the asymptotic solution of the stochastic equations ([K]), (|31|). 

Let us summarize the preceding discussion. The dynamical effects of 
the quenched disorder have been evinced by examining the varying stochastic 
features of a single random variable, namely, the probability density of the 
particle's occurence at the origin. Having a negative mean bias, the time- 
asymptotic and disorder-averaged value of this quantity is proportional to the 
absolute value of the mean force. In the Sinai-like case, i.e. for the vanishing 
mean bias, we have given the exact asymptotic formula describing the decay 
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of this quantity. The decay is slower than in the corresponding model without 
the reflecting boundary at the origin. Finally, having a positive mean bias, 
the particle escapes towards infinity, with a finite velocity or not, depending 
on the value of the trap-permeability parameter. The existence of deep traps 
with a long trapping time is crucial for the existence of anomalous dynamical 
phases. 

In the present work, we have chosen to formulate the diffusion problem in 
the presence of reflecting boundary conditions. After a slight modification, the 
method can be adapted to other types of boundary conditions. For instance, 
taking a fixed probability density at two boundaries, our method can be used to 
the analysis of the stationary-flux distribution in the one-dimensional random 
medium Jl8| . 



We have not aimed at the exhaustive description of the particle dynam- 
ics. Instead, we have concentrated on features which can be directly related to 
the probability density at the origin. The detailed description of the "noise" 
P(A; z) allows, at least in principle, for an investigation of other aspects, such 
as the time-asymptotic thermally-averaged first moment of the particle's po- 
sition. Another example would be the disorder-averaged diffusion coefficient. 
Its analysis requires the calculation of the second thermally-averaged moment 
of the particle's position on the one hand, and the higher-order (generally non- 
self-averaging) terms in the small- z expansion of the first moment on the other 
hand. We have shown that the thermally-averaged moments obey a system of 
stochastic differential equations with P(A; z) playing the role of the "input" 
noise. The probabilistic description of the moments will be reported elsewhere. 

Summing up, the paper presents an approximation-free study of the dif- 
fusive dynamics in an one-dimensional Markovian Poisson random potential. 
It provides a firm basis for the intuitive understanding of diffusion in more 
involved circumstances. 
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Appendix 



Let us consider the integrals 



h(z) 



clef 



3-W 

<2+0) 
g_0) 



g - 





(3) 


-q][q- 


- „_(*) 








- q+( z ) 




.5 





v+(z) 



[q - q+(z)][q - q+( z )] 



(A.1) 



Hence Iq(z) is equal to the integration constant C(z) in Eq. (|43|) . Having 
known /ol 2 ) and /^(z) for some k > 1, the fc-th stationary moment (Q k (z)) is 
simply given by the ratio Ik{z) / Iq(z). In this Appendix, we want to display 
some intermediate steps in the calculation of the above integrals. In order to 
keep the following formulae in a reasonable shape, we shall write v± instead of 
the more descriptive designation v±(z), used in the main text. Similar remarks 
hold for the four ^-dependent quantities q± and q'±. 
First, we notice the equality 



(g- -q)(q- qt) + (g - g+Xg - q'+) 

/+-/- 



(A.2) 



and we introduce the substitution x = (g — g + )/(g_ — q + ). Thereupon, the 
integrals Ik{z) assume the form 



h{z) 



dx [(1 — x)(x + a) + x(x + 6)] fc+1 x 



(/+-/-)* 
(1-x) 



u-—l x u+ ~ 1 



(x + a 



x + 6) 



(A.3) 



where a = (q + - q'_)/ (q_ - q + ) and b=(q+- q' + )/ (g_ - g+). 

Then, using the binomial theorem, we expand the (k + l)-th power in 
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(|A.3|) and write Ik{z) as a sum of k + 2 integrals: 



(u-f-) k UM k+1 -fi 

1 



dx (1 — x 



x + a 



Finally, the integrals Jk,j(z) are expressed through the Appell function T\ p5 
(the hypergeometric function of two variables |23|, 23). We have 



■4,j0) 



„ ++fc ,r( I /_ + j)r(i/+ + A;-j + i) 



r(i/ + A; + 1) 



X 



x (a, /?,/?', 7; 



-a , - 



(A.6) 



with a = z/ + + — j + 1, /3 = V- — j + 1, P' = u + — k + j, 7 = v + k + 1, 
and v = V- + u + . In the small- £ limit, the behaviour of the result depends 
on the two variables a and 6; one can employ the functional relations between 
the Appell functions. Moreover, having a special relationship between the four 
parameters, the Appell function can be expressed in terms of the ordinary 
hypergeometric function. In our context, the case k = is exceptional due 
to the possibility of such reduction. Alternatively, this point can be seen by 
introducing the second substitution, y = [x(l + b)]/(x + b). One then obtains 



Jo,o( z ) 



a u - +1 (l + b) u + +l Jo 

(1 + a)- 1 r(i/_)r(i/+ + 1 



a"- (1 + 6)"+ r(i/ + l 
1 n 



dyy'+il-yy-^l + uyr"- 1 
F(v~, v+, v + 1; —u) , 



(A.7) 



6a"- (1 + 6)"+ Jo 



dyy v +- 1 (l-y) v -(l + uyY 



6- 1 I>_ + l)r(i/+) 

— .7 r (z/_,z/ + ,z/+ 1; 



a"- (1 + 6)"+ T(z/ + 



(A.8) 
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with v = V- + u + and u = {b — a)/[a(l + b)]. Summing these two expressions, 
we acquire 



Io(z) 



Jo,o( z ) + J oA z ) 1 1 



q- - ^ 



Vt a v -(l + b)"+ 



B(v-, v + ) T{y-, v+, u + 1; -u) 



(A.9) 



where B(x, y) is the Euler beta function |23| and the factor 



„ n+(<7- - q'_) + n_(g + - q' + ) n + ^Az + g + n^Az + ft 



n_ + n + 



n_ + n + 



can be regarded as the mean value ^4z + <fi 2 (\)j — c.f. Eq. (f2~6|). Finally, 
upon inserting the definitions of the ^-dependent quantities a = a(z), b = b(z) 



and u = u(z), one recovers the normalization constant (45) from the main 
text. 
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Figure 1: 
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Figure 3: 



Force f(x) 

U- — 




Figure 1: Typical realization of the Markovian Poisson dichotomic force 
taking two values of opposite signs, /_ < < /+ with / + = 2|/_|. The slope 
of the corresponding potential is alternately positive and negative: there are 
traps in this case. 




Figure 2: Typical realization of the Markovian Poisson dichotomic force taking 
two values of the negative signs, /_</+< with |/_| = 2\f + \. The slope of 
the corresponding potential is always positive. There are no traps in this case, 
the particle being stuck towards the reflecting boundary at the origin. 




Figure 3: Typical realization of the white shot-noise random force and of the 
corresponding potential. The lengths of the vertical segments of the potential 
are independent, exponentially distributed random variables, their mean value 
being 7. There are traps in this case. 



